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I. INTRODUCTION 


For a fluid supercooled isobarically below the melting temperature T m , the ordered (say, 
crystalline) state is thermodynamically favorable. At the moderate supercooling (T m — 
T)/T m , the transition into an ordered state is started through nucleation mechanism, that 
involves emergence of the crystalline nuclei, which are able to grow; T is temperature. Hence, 
the behavior of the overall transition should be essentially determined by the rate charac¬ 
teristics: the waiting time of the first crystalline critically-sized nucleus Ti, the nucleation 
rate J that is amount of the supercritical nuclei formed per unit time per unit volume, and 
the growth rate v g that specifies growth law of the supercritical nuclei. 

Some general features of the nucleation kinetics can be comprehended within the classical 
nucleation theory- - - and its extensions (see in Refs.- - -). According to the classical view, 
the driving force for the nucleation grows upon the increase of supercooling. This means 
that with the increase of supercooling an ordered state becomes thermodynamically more 
favorable, while the waiting time for nucleation T\ and the nucleation time scale 1 / J must be 
shortened. On the other hand, with temperature lowering (below the melting temperature 
T m ), the mobility of molecules (atoms) decreases. As a result, any structural rearrangements, 
including these responsible for nucleation, must be suppressed by the growing viscosity. 
When a fluid is cooling down without crystallization to temperatures corresponding to the 
viscosity r](T g ) > 10 12 - 10 13 Pa-s, it is “freezing” as disordered solid; where the temperature 
T g is identified with the glass transition temperature. Although crystallization of glasses 
proceeds over time-scales^, which are commonly larger than experimentally acceptable, 
the structural ordering in a glass can be accelerated by out-of-equilibrium processes resulted 
from reheating or applied shear deformation^ - —. On the other hand, there are indications 
(see Refs.— - —) that the time-scales of structural relaxation and of ordering in glassy systems 
becomes shorter, when we move over equilibrium phase diagram to the range of more higher 
pressures. 

Moreover, debated issues in the field are related to the temperature dependence of the 
transition rate characteristics at deep levels of supercooling^ - —. So, for example, empir¬ 
ical T-dependencies of the nucleation lag-time and of the steady-state nucleation rate are 
discussed in review^, where results for some stoichiometric glasses (3Mg0-Al 2 0 3 • 3Si0 2 , 
Li 2 0-2Si0 2 , Na 2 0-2Ca0-3Si0 2 , et al .) are given. As it is demonstrated in Ref. 2 - within the 
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available experimental data, the lag-time of nucleation and the steady-state nucleation time- 
scale 1 / J s reach the lowest values at a certain moderate levels of supercooling. Moreover, 
both rate terms start to grow with the further increase of supercooling and with the ap¬ 
proaching the glass transition temperature. Remarkably, the possible correlation discussed 
in Ref.— between some features in the temperature dependencies of these rates (for example, 
the maximum steady-state nucleation rate) and the reduced temperature T g /T m provides, 
in fact, indirect implications about “unified laws”, which can be inherent in the nucleation 
kinetics. In this work, we extend this view by focusing on the crystal nucleation time T\, 
identified here as the average waiting time for the first critically-sized nucleus^ 1 ^. 

The possibility of unified description using scaling relations has been proposed and stud¬ 
ied for the case of nucleation of liquid droplets in the condensation process. Here, an in¬ 
triguing feature emergent in the analysis of data for the vapor-to-liquid nucleation is a 
supersaturation-temperature scaling of the nucleation rate data^—. Namely, as shown by 
Halevi - — , data for the nucleation rates plotted vs. C' 0 In S/[T C /T — l] 3 / 2 can collapse onto 
a single line. Here, S = p/p C oex is the supersaturation, p is the pressure of supersatu¬ 
rated vapor, Pcoex is the pressure at the coexistence curve, T c is the critical temperature, 
and Cq is a normalization factor. This result is very interesting for the following reasons. 
First, scaling relation allows one to compare the nucleation data of various independent 
studies for a system, even though those studies have not the identical pressure-temperature 
(supersaturation-temperature) conditions. Moreover, if the scaling is valid, then this is 
indication that there is a single reduced variable instead of the pair, T and S'; and this 
variable is sufficient for a unified description of the steady-state vapor-to-liquid nucleation 
rate. Recently, Diemand et al.—— suggested a new scaling relation for the nucleation rate of 
homogeneous droplets from supersaturated vapor phase, where other set of the parameters 
is utilized. From the mentioned considerations, it is reasonable to try to extend the ideas of 
scaling relations to the case of other transition — to the case of crystallization^. The present 
study is mainly aimed at the consideration of this issue. For this, the analysis of the crystal 
nucleation times from the experiments and molecular dynamics simulations is carried out 
for several systems at temperatures T <T g . 

The paper is organized as follows. In Sec. [Hi the reduced temperature scale is introduced. 
Section m presents the simulation details and computational methods. It includes the 
description of two model systems taken for molecular dynamics simulations, the cluster 
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analysis and the statistical method utilized for the evaluation of the nucleation characteristics 
within the simulation data. Discussion of the results is given in Sec. HVI The main conclusions 
are finally summarized in Sec. |V] 


II. REDUCED TEMPERATURE SCALE 

In evaluation of the unified temperature dependencies for characteristics of the super¬ 
cooled liquids, one encounters the problem that the interested temperature range, 0 < T < 
T m , contains three control points - the zeroth temperature T = 0 K; the glass transition 
temperature T g and the melting temperature T m , - where T g and T m in the Kelvin scale 
have not the same values for different systems. Therefore, there is necessity to use a reduced 
temperature defined usually either through T m , or through T g , depending on the problem 24 . 

For example, according to Angel l' 37 " 38 , the inverse reduced temperature T g /T is used to 
fulfill the “strong-fragile” classification of viscous (supercooled) liquids by means of the plot, 
in which the viscosity in logarithmic scale, log 77 , is considered as a function of T g /T. Since 
one has log[? 7 (T g )] = 12-13 for all the supercooled liquids by definition, then the values of 
log 77 will be comparable on the reduced temperature scale 0 <T g /T < 1 in neighborhood 
of T g . Further, the supercooling (T m — T)/T m or its conjugate quantity T/T m represent 
also the reduced temperature scales (Ref.—) and are used to compare the characteristics 
of supercooled liquids for temperature range T < T m . Here, a reasonable consistency is 
ensured as we approach the melting temperature T m . Ambiguity of the choice of reduce 
temperature scale is because the ratio T g /T m depends on the system (material) and can 
be different even for the systems of same type. For example, the ratio of T g /T m for glasses 
Li 2 0-2Si0 2 , Ba0-2Si0 2 and 2Na 2 0-Ca0-3Si0 2 , which belong to the group of silicate glasses, 
does not have the same value, and is equal to 0.56, 0.568, 0.512, respectively—. Moreover, 
the quantity T g /T m is dependent on cooling rate dT/dtt applied to prepare glass at a desirable 
temperature and can have different values for the different isobaric lines of a phase diagram. 
Therefore, the absolute temperature T as well as the reduced temperatures T/T g and T/T m 
can not be considered as convenient parameters, with respect to which evaluation of the 
unified regularities could be examined. 

To overcome this one needs to specify a temperature scale T, in which the control points 
mentioned above - the zeroth temperature, the glass transition temperature and the melting 
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temperature - are fixed and have same values for all systems. We suggest a possible simple 
way to realize this. Let us define the following correspondence between the values of T for 
the three temperatures (the zeroth temperature T = OK, the glass transition temperature 
T g and the melting temperature T m ): 


T = 0 

Tg = 0.5 
T m = 1 


at T — 0 AT, 

(la) 

at T — T g , 

(lb) 

at T — T m . 

(lc) 


The conditions (PQ) are fulfilled with the simple parabolic relation: 
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With the known T m and T g for a system, relation (|2]) provides transform of the absolute 
temperature scale T into the reduced scale T, where all the temperature points coincide for 
all considered systems (see Fig. [□). Moreover, when the ratio T g /T m approaches value 0.5, 
the quadratic contribution in Eq. ([2]) vanishes and Eq. (J2J) is simplified to 

1 T 


T ~ 


2 T n 


It should be pointed out that the temperatures T = 0 K, T g and T m in relation ([2]) 
correspond to the same isobar. Relation (J2]) transforms the (p, T ) phase diagram within 
the range 0 < T < T m to the (p, T ) phase diagram unified for all systems, where the glass 
transition line and the melting line are parallel to the ordinate, p-axis, and intersect the 
abscissa at T g = 0.5 and T m = 1, respectively. 
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FIG. 1. Demonstration of the transformation of the absolute temperature scale T into the reduced 
temperature scale T, which is system-independent and characterized by fixed values of the melting 
temperature T m = 1 and the glass-transition temperature T g = 0.5. Simple realization of the 
transition can be done by means of relation 

III. SIMULATION DETAILS AND COMPUTATIONAL METHODS 

In this work, we consider two model systems - the Dzugutov (Dz) system^' 40 . an d the bi¬ 
nary Lennard-Jones (bLJ) mixture^—. Both systems are known as the model glass-formers 
suitable to study the properties of glasses by means of molecular dynamics simulations^ - —. 
In this work, the glassy samples were generated by fast quench of equilibrated fluids at the 
fixed pressure p. The corresponding pathways are shown on the phase diagrams in Fig. [2l 
Consideration of the (p, T)-points of the phase diagrams with high values of the pressure p 
allows us to deal with such conditions at which the structural ordering in the glassy systems 
proceeds over time-scales available for simulations even at temperatures below T g (Refs.- 1 ^). 
Thereby, the value of the pressure p was chosen so that a clearly-detected nucleation event 
was observable over the simulation time scale. Hence, the simulations with the generated 
glassy samples were performed in the ApT-ensemble; N is the number of particles. The con¬ 
stant temperature and pressure conditions are ensured by using the Nose-Hoover thermostat 
and barostat. For each (p, T)-point, more than fifty independent samples were generated, the 
data of which were used in a statistical treatment. For a single simulation run, N = 6 912 
particles were enclosed in a cubic cell with periodic boundary conditions. Note that the 
terms e and a define the units of energy and length, respectively. Time, pressure, and 
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temperature units are measured in ro = a\Jm/e , e/cr 3 , and e/ks, respectively. 




FIG. 2. Pressure-temperature phase diagram for the Dz-system (left panel) and for the bLJ- 
system (right panel). The full curves denote the boundary between liquid and solid phases; the 
curve for the Dzugutov system is reproduced from data of Fig. 4 in Ref. 13 . The dashed curves 
mark the boundary between the supercooled liquid and the amorphous solid, when liquid is cooled 
during isobaric simulations with the rate dT/dtt = 0.001 £/(/cbTo). The full squares indicate the 
equilibrium liquid states, which were used as starting points to generate glassy samples. Pathways 
related with preparation of the glassy samples are schematically shown by dotted arrows; and the 
full circles denote the (p,T)- points, at which the transition into ordered states was tracked. 


The Dzugutov system. - In case of the Dzugutov system, all particles are identical and 
interacting via a short-ranged pair potential 


U Dz ( 


= A(r* 5 
+ B exp 


— B ) exp 
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0(a — r*) 
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r* = r AL. 
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(4) 


where @(...) is the Heaviside step function, the values of parameters A = 5.82, B = 1.28, 
m — 16, a = 1.87, 6 = 1.94, c = 1.1 are chosen as suggested originally in Ref.—. The 
simulations were performed for the system along the isobaric line with the pressure p = 
14 e/a 3 at the temperatures T = 0.05, 0.1, 0.15, 0.3 and 0.5 e/ks below T g ~ 0.65 e/ks- For 
the isobar, the melting temperature is T m ~ 1.51 e/ks, that yields the temperature ratio 
T g /T m ~ 0.43. 
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The binary Lennard- Jones mixture. - The semi-empirical (incomplete) Lorentz-Berthelot 
mixing rules^k—, 

&bb = 0.8 a aa-, 

O'AA + &BB 

° AB - 2 -’ 

£bb = 0.5 £aa, 

£ab = £aa + £bb, 

were utilized at the simulations of the binary Lennard-Jones system A^qB^q with the poten¬ 
tial 
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where a, /3 E { A , B}, the labels A and B denote the type of particles, r ig is the distance 
between the centers of particles i and j. Note that we take e = £aa, & = &aa, and the 
mass of a particle is m = tua = mg = 1. For the bLJ system, we consider the isobar with 
the pressure p = 17e/a 3 at the temperatures T = 0.01, 0.05, 0.1, 0.2 and 0.3 e/ks, that 
are lower than the transition temperature T g ~ 0.92 e/fcg. The isobar contains the melting 
point with T m ~ 1.65 e/k,B- Therefore, the temperature ratio is estimated as T g /T m ~ 0.56. 

Cluster analysis. - The local domains of a crystalline symmetry are examined by means of 
the cluster analysis^ 1 ^, introduced originally by Wolde-FrenkeP 9 . The consideration of the 
local environment around each particle is performed by means of 13-dimensional complex 
vector with the components^ 


N b (i ) 


%r, 


U = 


•\ 'y 1 Pij). 


( 6 ) 


Here, Ye >rn {dij-,Tij) are spherical harmonics, lV&(i) is the number of neighbors for i particle, 
dij and ipij are polar and azimuthal angles, which characterize the radius-vector r t j. Then, 
the local order for each i particle can be numerically evaluated by means of the parameter— 

(a 6 \ 1/2 

= bj E IfcMI 2 • (?) 

\ m =—6 / 

whereas degree of the orientational order can be estimated by means of the global orienta¬ 
tional order parameter Q e defined as an average of q 6 (i) over all N particles^: 
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For a fully disordered system the parameter q%{i) is close to zero, while it grows with increas¬ 
ing structural ordering. For perfect fee, bcc and hep systems one has the largest possible 
values for the parameters^: 


9e(*) 

= Q e ~ 0.5745 

(fee), 

Qe(i) 

= Qe ~ 0.5106 

(bcc), 

Qe(i) 

= Q 6 ~ 0.4848 

(hep). 


First, we define “neighbors” as all particles located within the first coordination, the 
radius of which is associated with position of the first minimum in the pair distribution 
function^. Further, according to the Wolde-Frenkel scheme^ we specify the pair of neigh¬ 
boring particles (i and j) as connected by a crystal-like bond if the following condition is 
fulfilled: 
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Condition (J9J) allows one to distinguish the particles correlated into an ordered structure^. 
Finally, particle i is identified as included into a crystalline structure if it has four and 
more crystal-like bonds. The last condition is applied to exclude from consideration the 
structures with a negligible number of bonds per particle, which occurs even in equilibrium 
liquid phased By means of this routine, the particles involved into the crystalline domains 
are detected. 

Figure [3] demonstrates, as an example, the crystalline clusters emerging in the glassy 
Dz-system at T — 0.5 e/ks over the transient nucleation regime, where no nuclei capable 
to grow are detected [Fig. [31(a)], and at the time t = 250 r 0 , when the first nucleus of the 
critical size appears [Fig. [31(b)] . 


Statistical treatment of the cluster analysis results. - The growth trajectories of the crys¬ 
talline nuclei, n a i(t), extracted from the different simulation runs are treated within the 
mean-first-passage-time method^—. Here, n defines number of the particles involved in the 
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FIG. 3. (Color online) Snapshots of the Dz-system at T = 0.5 e/ks and at different times, for 
which the particles recognized as belonging to the crystalline phase are shown only, (a) System at 
the transient nucleation period; t = 100 to- There are no nuclei capable to grow, and their sizes are 
smaller than the critical size n c . (b) System at the time t = 250 To, when the first critically-sized 
nucleus emerges; n c — 105 partilces. The critically-sized nucleus is marked by red circle. 



FIG. 4. (Color online) Mean-first-passage-time distributions ri(n) and its first derivatives 
dri(n)/dn defined from simulation data for the Dz-system and the bLJ-system at different tem¬ 
peratures. 
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nucleus at the time t, the mark a denotes the index of simulation run, whereas the order 
number of the nucleation event i indicates that the it\i nucleus of the size n appears at the 
time t during the ath simulation run. On the basis of the extracted trajectories n a i(t), the 
mean-first-passage-time distributions Tj(n) are evaluated for each ith-order nucleus (for de¬ 
tails, see Ref.—). Further, the critical size n c and the average waiting time for the All-order 
nucleus, r tl i — 1, 2, ..., are defined from the analysis of the distributions Tj(ra) and of the 
hrst derivatives dTi(ri)/dn, according to the scheme suggested in Ref.—. In this work, we 
focus on the characteristics for the largest nucleus — i.e. on its critical size n c and average 
waiting time T\. 

As an example, we show in Fig. 0] the mean-first-passage-time distribution Ti(n) and its 
hrst derivative dri(n)/dn computed for both the systems. As can be seen, the distributions 
Ti(n) are characterized by three regimes. The hrst regime, for which small values of n 
correspond to Ti(ri) with zero value, is associated with pre-nucleation. Here, the nuclei with 
different sizes (albeit, small sizes) appear with equal probability. The second regime, in 
which the distribution T\ (n) has the pronounced non-zero slope, contains information about 
a nucleation event. Namely, detected from the hrst derivative dr\{n)/dn location of an 
inflection point in the distribution T\ (n) for the regime defines the critical size n c , whereas 
Ti(n c ) = T\ is directly associated with the average waiting time of the hrst critically-sized 
nucleus^. Finally, the third regime, where the slope of T\{n) decreases, corresponds to 
growth of the nucleus. Note that such shape of the mean-hrst-passage-time distribution is 
typical for an activated process. The absence of the pronounced plateau in Ti(n) for the 
third regime indicates that the nuclei growth proceeds over a time-scale comparable the 
nucleation time T\ (Ref.—). 


IV. DISCUSSION OF RESULTS 

We start from evaluation of some properties of the nascent ordered structures, that can 
help to elucidate the mechanism of the ordering. Figure 0 shows the time-dependent order 
parameters - the global orientational order parameter, and the size of the largest 

cluster, n(t), - evaluated on the basis of the simulation data. In initial stages, the param¬ 
eters Q${t) and n{t) fluctuate around their starting values. After an incubation time, both 
the parameters start to growth rapidly. Such evolution of the order parameters indicates 


11 



FIG. 5. Trajectories of the global orientational order parameter Qo(t) and of the largest crystalline 
nucleus size n(t) defined from a single simulation run for the Dz-system (left panel) and the bLJ- 
system (right panel). The dotted horizonal lines on the plots for n(t) correspond to the critical 
sizes n c defined from the statical analysis within the mean-first-passage-time method. 


on activated character of the transition^. The nucleation event is clear detectable on a 
particular trajectory n(t), where it is associated with the start of sharp grow of n(t). While 
rough estimates for the nucleation time-scale T\ and for the critical size n c can be done even 
from the particular trajectories n(t) [see Fig. [5], the averaged values for both the quantities 
can be computed directly by means of the statistical method presented in Sec. IIII1 

Further, cluster analysis reveals that the nuclei of the critical size are localized. As 
contrasted to Ref.—, no ramified structures were detected even at very deep supercooling. 
For quantitative characterization, the asphericity parameter So was computed according to 

(j - I ) 2 + (I - I ) 2 + (I -I ) 2 ' 

y ± xx ± yy) ' y ± xx ± zz) ' y ± yy ± zz) 


s n = 
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( 12 ) 


defines the components of the moment of inertia tensor associated with a critically-sized 
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nucleus; the brackets (...) mean the statistical average over results of the different simulation 
runs. The parameter So approximates the unity, So —> 1, for an elongated and ramified 
cluster, and one has So —> 0 for a cluster, the envelope of which is of spherical shape. For 
both the systems (Dz and bLJ), we find that the asphericity parameter is So — 10 3 , and 
the size of the critical nucleus remains finite. This is evidence that the transition into an 
ordered phase is initiated rather through nucleation mechanism, and that is in agreement 
with findings of Refs.—' For the Dz-system, our estimations reveal that the critical size 
changes from n c = 108 ± 5 to 88 ± 6 particles with the temperature decrease (increase of 
supercooling) from T = 0.5 e/kB to 0.05 e/fcg. For the bLJ-system, we find that the critical 
size decreases from n c = 59 ± 4 to 42 ± 3 particles with the temperature decrease within the 
range 0.3 e/k B >T> 0.01 e/k B . 

Figure [6] shows the values of the average waiting time of the first nucleus of the criti¬ 
cal size, Ti, estimated from simulation data for the Dz-system and the bLJ-system at the 
different temperatures. We note that the deep levels of supercooling are considered for 
both the systems corresponding to the temperatures much below T g . The particle mobility 
diminishing with supercooling results in the growth of T\ with the temperature decrease. 
The finite values of T\ comparable with the duration of numerical experiment may seem 
surprising for a glassy system. Actually, the microscopic kinetics of a glass changes with 
moving over phase diagram for the range of high pressures^. Namely, at high pressures the 
structural relaxation as well as the transition of glassy system into a state with the lower 
free energy proceeds over shorter time scales^——. Therefore, the reduction of the values of 
t i is admissible for the range of phase diagrams. 

Although the quantity T\ for both the systems demonstrates similar temperature de¬ 
pendence, it is difficult to say something about quantitative correspondence to the general 
nucleation trends. Is such temperature dependence of the nucleation waiting time, ri(T), is 
typical for the considered thermodynamic range or not? One of the possible ways to clarify 
this is to bring the extracted values of the nucleation waiting time t\(T) into a unified scaled 
dependence. To construct scaling relation, we propose to use the reduced temperature T 
defined by relation (|2|) , in which the values of the glass transition temperature and the melt¬ 
ing temperature are fixed for all systems. Then, the simplest nonlinear T-dependence of iq 
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FIG. 6. (Color online) Average waiting time of the first critically-sized nucleus t\ (in units of to) 
versus reduced temperature for the Dz-system ( T g /T m = 0.43) and for the bLJ-system ( T g /T m = 
0.56). The spanned thermodynamic ranges correspond to deep levels of supercooling with the 
temperatures below T g . 


can be chosen in the form: 



(13) 


where rf is the average waiting time for the first critically-sized nucleus at the state with the 
temperature T g (we remind that T g = 0.5). The dimensionless parameter 7 > 0 characterizes 
ability of the system at the considered (p, T)-state to retain structural disorder. In particular, 
the exponent 7 takes high values for the system with good glass-forming properties, and must 
be characterized by small values for the fast crystallizing systems. Since the nucleation 
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TABLE I. The melting temperature T m , the ratio T g /T m , the waiting time for the first critically- 
sized nucleus rf at the transition temperature T g , the exponent 7 estimated from Eq. (fl3l) . the 
parameters K\ and K 2 evaluated by Eq. (I3bl) for several systems. 


System 

T m 

Tg/Tm 

rf 

7 

K\ 

K, 

Dz (at p = 14e/cr 3 ) 

1.51 s/ks 

0.43 

211 TO 

0.27 

0.553 

-0.053 

bLJ (at p = 17e/a 3 ) 

1.65 e/ks 

0.56 

760 t 0 

0.025 

0.427 

0.073 

Li 2 0-2Si0 2 

1286 K 

0.56 a 

1869 seel 

70 

0.424 

0.076 

Na 2 0-2Ca0-3Si0 2 

1549 K 

0.53^ 

5150 seel 

50 

0.466 

0.034 

K 2 0-Ti0 2 • 3Ge0 2 

1308 K 

0.631 

990 sec £ 

30 

0.281 

0.219 


a Experimental data of Ref.—. 
b From experimental data of Ref—. 
c Experimental data of Ref—. 
d From experimental data of Ref—. 
e Experimental data of Ref—. 
f From experimental data of Ref—. 


waiting time 17 varies with pressure, then the exponent 7 should be dependent on the 
pressure, at which a supercooled liquid evolves. Namely, the exponent 7 is decreasing 
function of the pressure p for the systems, in which the nucleation time scale decreases with 
pressure. The Dz and bLJ systems correspond to the case. 

To verify validity of relation (fTBjh we place the rescaled data for the average waiting time 
t 1 for the Dz and bLJ systems vs. the reduced temperature T on the common Fig. [0 For 
clarity, the axis of ordinates is presented on a logarithmic scale, where the fitting parameter 
7 corrects the slope in accordance with the master-curve 


Tl 


7 



(14) 


which appears from (fTcfll at the exponent 7 = 1 . The reduced temperature T in Eq. 
guarantees that the temperature points spread over the abscissa in the same manner for 
all the considered systems, whereas the dimensionless parameter 7 forces all the ordinate 
points to collapse onto the master-curve (fT4)h Since our simulation results for the Dz and 
bLJ systems cover the temperature range T < T g and we did not estimate the nucleation 
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time t i at the transition temperature T g , then the term rf was taken as a fitting parameter. 
Namely, its values were found by extrapolation of the data for T\ to the temperature point 
T g = 0.5, where the function (l/y) log(ri/rf) must be equal to zero (see Fig. [7]). Numerical 
values of rf are given in Tab. |U As can be seen from Fig. [TJ all the data obtained on the basis 
of molecular dynamics simulations follow the unified master-curve. Moreover, in contrast to 
the case of the Kelvin temperature scale, values of the melting temperature T m and of the 
transition temperature T g are not dependent on pressure. Therefore, the results shown on 
Fig. [7] can be supplemented by the data for any supercooled liquid at arbitrary value of the 
pressure p. 

Moreover, it is attractive to extend the study and to verify the scaling law (IT3]) with 
the experimental data. While the direct experimental measurements of 77 are difficult^, 
we suggest the next routine for the approximative estimation of 77, which can be realized 
with the experimentally measurable quantities - the steady-state nucleation rate J s and the 
induction time r m ^. According to Kashchiev-— , the number density of the supercritical 
nuclei in the system, iy, evolves with time as 

iy(t) t 


TstJind 7ind 


- 1 


(15) 


llA(-l)' 

7T 2 
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77r 


exp 


m 2 n 2 t 


6 Ti. 


ind 


For the time t — 77 one has iy(ri) = 1/K, and Eq. (fl5|) takes the form: 
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JstV 
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Further, Eq. (ITTH) was numerically solved with the experimental J s and T; n d for Li 2 0-2Si02 
reported in Ref.—, for Na 2 0-2Ca0-3Si0 2 presented in Ref.—, and for K 2 0-Ti0 2 -3Ge0 2 given 
in Ref.—. The extracted rescaled values of the average waiting time 77 are also presented in 
Fig. [0 As can be seen from Fig. [7] “experimental” data for t\ provide the T-dependence, 
which is in agreement with scaling relation (I13p as well as with the simulation results for 
the bLJ-system and the Dz-system. 

Analysis of the reduced temperature scale T for the systems reveals that the quadratic 
contribution in equation for T [see Eq. (|2])] can be insignificant as for the Dz-system and for 
Na 2 0-2Ca0-3Si0 2 , where the ratio T g /T m is equal to 0.43 and 0.53, respectively (see Tab.[J). 
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FIG. 7. (Color online) (a) Scaled waiting time for the first critically-sized nucleus ( 1 / 7 ) log 10 (ri/rf) 
is plotted as a function of the reduced temperature T. Here, T g = 0.5 is the scaled transition 
temperature (marked by arrow), rf is the waiting time at the transition temperature T g , and 7 is 
the htting parameter, (b) The same but for the temperature range 0.46 <T< 0.58. Values of the 
parameters rf and 7 are given in Tab. [H Consistency of the data to the master-curve, which set 
the temperature-dependence T g /T, provides support to the validity of scaling relation (fT3l) . 


With away from value 0.5 for the ratio T g /T m , weight of the quadratic contribution, K 2 , 
increases. The values of the parameters K\ and K 2 are comparable for K 20 -Ti 02 • 3 GeC >2 
characterized by the ratio T g /T m = 0.63. To our knowledge, the highest value of the ratio 
T g /T m appears for Na 20 -Al 2 03 • 6 SiC >2 and is equal to 0.78 (Ref.—). 

The values of the exponent 7 differ for the considered systems by four orders of magnitude, 
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and an order of magnitude between the Dz-system and the bLJ-system (see Tab. [[])• The 
large scatter in the values of 7 is due to the change of the waiting nucleation time 77 within 
the temperature range differs essentially for the systems. One can demonstrate this with 
the results for the Dz and bLJ systems shown on Fig. O Within the temperature range 
0.025 < T/T m < 0.2 the time scale 77 is changed by the factor 0.625 for the Dz-system, 
whereas it changes by the factor 0.96 for the bLJ-system. For the systems with complicated 
structural units (i.e. for silicate glasses) the change is much more pronounced^. The case 
with the smaller change in the temperature dependence of 77 will corresponds to the smaller 
values of the exponent 7 in scaling relation (f]ySjh 

Although the scaling law HTTP is suggested rather as an empirical result, its qualitative 
justification can also be done. At the temperatures comparable with and lower than the glass 
transition temperature T g , the local structural rearrangements responsible for the nucleation 
are driven rather by kinetic aspects associated with the viscosity than by thermodynamic 
contribution. Therefore, it is reasonable for the range of high supercooling to expect the 
existence of correlation between the waiting time for nucleation 77 (F) and the structural 
relaxation time r a {T ) ~ 77 (F), and, thereby, between the time 77 (T) and the viscosity 77 (F): 

ri (T)~r,(T). (17) 


Hence, the Vogel-Fulcher-Tammann equation provides the most popular viscosity model 
(this equation is also known as the Williams-Landel-Ferry rnodel^—): 

l«g 10 V(T) = logm Voc + T , (18) 

where T c is the critical temperature of this model. Another equation for viscosity similar to 
VFT-model is provided by the mode-coupling theory——: 


V(T) = 


''lc 


7m 


(19) 


(F — Tmct) 

where Tmct is the (critical) mode-coupling temperature. The parameters r /oo, A and 
take positive values and are obtained by fitting Eqs. (fT 8 |h (HUD to experimentally measured 
viscosity data£^. Both the models predict a divergence of the viscosity 77 (F) when F — y T c 
(and F —> T M ct)- Moreover, both the models are able to reproduce 77 (F) for the supercooled 
liquid phase, i.e. T > T c (and F > Tmct), and are not applicable for the temperature range 
below F c (below Tmct ) because of a divergence in the temperature dependencies. On the 
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other hand, for a high-viscosity regime corresponding to the temperatures T < T g , the 
experimentally measured temperature-dependence of the viscosity rj(T) is reproducible by 
the Arrhenius law (see, for example, Fig. 6 in Ref.—), which is generalized by the Avramov- 
Milchev equationSA^: 


or 


l °SioV(T) = logio^ + 


logio 


'viry 

. Voc _ 




where A and a! are positive. 

On the other hand, let us now reconsider scaling relation (fl3j) . which can be rewritten in 
the form 

^=( 2 fy \ ( 22 ) 

since T g = 0.5. After substitution of Eq. d2j) into relation (l22l) and using the expansion 
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we obtain for the temperature range 0 < T < T m the following equation: 
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Assuming that proportionality in Eq. (1T71) 
. A simple analysis reveals that Eq. (1241) is 


where the parameter K\ is defined by Eq. 
holds, one can compare r.li.s. of Eqs. (l2Tj) and 
able to approximate the power-law dependence of Eq. (l2Tj) and generalizes the temperature 
dependence for the viscosity given by the Avramov-Milchev equation. Thereby, the scaling 
relation (IT3]) and the viscosity model with Eq. (12T| can be considered as consistent. 

Moreover, the fragility of a system can be estimated by means of the index m defined 


as^ 


m = 


d logio (v) 


T=T„ 


8{T g /T) 

Then, from Eqs. (ITT)) . Q22j) and (125)) we obtain the following relation 

<9 lo gio( r i) 
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d(T,/T) 
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which after substitution of Eq. (j3b|) can be rewritten as 


ffere the contribution in square brackets is positive for the range 0 < ( T g /T m ) < 1. Last 
two relations indicates that the exponent 7 and the index m are correlated terms, whereas 
7 can provide an estimate of fragility. 

V. CONCLUSION 

The mechanism of the structural ordering in the supercooled melts at extremely deep level 
of supercooling is one of the most debated issues in the consideration of the crystallization 
kinetics^ - - 0 . Let us mention some viewpoints in this regard. The mean-field theories, 
starting from the gradient theory of Cahn-ffilliard, provide indications that the structural 
ordering at a deep level of metastability can proceed through the spinodal decomposition— 1 . 
Interestingly, Trudu et al. for the freezing bulk Lennard-Jones system found a spatially 
diffuse and collective phenomenon of nucleation at deep supercooling. Authors treated such 
features as indirect signatures of a mean-field spinodal^. This was later criticized by Bartell 
and Wu^. According to experimental 72 and other simulations^ studies, the size of the 
critical embryo remains finite with decrease of the temperature of the supercooled liquid, 
in contrast to the mean-field theory predictions for a spinodal. Moreover, results of Ref .' 1 
reveal that crystallization in hard sphere glasses proceeds due to “a chaotic sequence of 
random micronucleation events, correlated in space by emergent dynamic heterogeneity”, 
and agree with findings of Bartell-WuS. In view of this, it remains still desirable to examine 
the mechanisms of the structural ordering in glasses within the new experimental/simulation 
results. 

In the present work, two model glassy systems with different interparticle interaction - the 
single-component Dzugutov system and the binary Lennard-Jones system - are simulated 
with the aim to study the structural ordering at deep supercooling. Remarkably, the simu¬ 
lation study covers a wide temperature range: from the temperatures comparable with T g to 
the temperatures corresponding to very deep levels of supercooling {T m — T)/T m ~ 0.97. By 
means of cluster analysis, we show that the structural ordering even at deep supercooling 
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proceeds through the formation of the localized crystalline domains, where the size of the 
critical embryo still remains finite. This supports the nucleation scenario of crystallization 
in the glassy systems, and is in agreement with the recent findings of Saika-Voivod et al. 
and Sanz et al.—. 


The average nucleation time is the quantity of main interest in the characterization of the 
initial stages in the nucleation kinetics. Here, it is estimated on the basis of the molecular 
dynamics simulation data (for the two model glassy systems) and from the available exper¬ 
imental data for the several glasses within the Kashchiev’s approximative equation. Our 
results show that, with the decrease of the temperature, the nucleation time T\ increases 
but still remains finite. Further, we find that the nucleation time T\ plotted as a function 
of the proposed reduced temperature follows the power-law dependence, unified for all the 
considered systems. The correlation between the proposed reduced temperature dependence 
for t i and the viscosity models for the amorphous solids supports the conclusion about the 
kinetic character of the initiation of the structural ordering in glasses, where the inherent 
glassy microscopic dynamics is predominating over thermodynamic aspects. 

Results of this study extend the idea of a unified description of the nucleation kinetics 
using scaling relations, which was originally applied to the analysis of the droplet nucle¬ 
ation rate data for the vapor-to-liquid transition (see^ and references to^). The later 
treatment indicates that the nucleation rate can be well described by the scaling function 
In {p/pcoex)/[T c /T - 1] 3/2 . In this study, we pursue a similar approach applied to crystal¬ 
lization and define such a variable, which might provide consistency in comparison of the 
crystal nucleation time data for different systems. Our realization differs from the scalings 
of Ref.—it is based on the reduced temperature scale with the fixed control points: the 
temperature T = 0, the glass transition temperature T g = 0.5 and the melting temperature 
T m = 1 for a considered system. Using this approach we find a correspondence of the scaled 
nucleation times as extracted from simulation and experimental data for the various systems 
to a unified power-law dependence. Finally, we note that because of experimental difficulties 
in extraction of the quantitative information about the initial stages of the crystallization ki¬ 
netics, few of the experimental studies cover the range of supercooling (T m — T)/T m > 0 . 6 ^. 
In this regard, it could be desirable to verify the suggested scaling law with additional 
experimental studies, especially, for the glassy systems at deep supercooling. 
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